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Investigations of the intrinsic electromagnetic nonlinearity of superconductors give insight into the 
fundamental physics of these materials. Phase-sensitive third-order harmonic voltage data its/ = 
\u3f\exp(i(j}3f) are acquired with a near-field microwave microscope on homogeneous YBa2Cu307_i 
(YBCO) thin fihns in a temperature range close to the critical temperature T^. As temperature 
is increased from below Tc, the harmonic magnitude exhibits a maximum, while the phase, 7r/2 in 
the superconducting state, goes through a minimum. It is found that samples with doping ranges 
from near optimal {5 = 0.16) to underdoped {S = 0.47) exhibit different behavior in terms of 
both the harmonic magnitude and phase. In optimally-doped samples, the harmonic magnitude 
reaches its maximum at a temperature Tm slightly lower than that associated with the minimum 
of phase Tm and drops into the noisefloor as soon as Tm is exceeded. In underdoped samples Tm is 
shifted toward lower temperatures with respect to r,„ and the harmonic voltage magnitude decreases 
slower with temperature than in the case of optimally-doped samples. A field-based analytical 
model of ua / is presented, where the nonlinear behavior is introduced as corrections to the low-field, 
linear-response complex conductivity. The model reproduces the low-temperature regime where 
the (T2 nonlinearity dominates, in agreement with published theoretical and experimental results. 
Additionally the model identifies Tm as the temperature where the order parameter relaxation 
time becomes comparable to the microwave probing period and reproduces semi-quantitatively the 
experimental data. 
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I. INTRODUCTION 

Recently, a number of experiments have shown evi- 
dence of unusual properties above the superconducting 
transition temperature in under-doped cuprate supercon- 
ductors. Observations include a significant Nernst effect 
in the pseudo-gap phase of La-Sr-Cu-0, suggesting the 
existence of vortex excitationsi. Diamagnetic response 
above Tc has also been observed in under-doped Bi-Sr- 
Ca-Cu-0^. These properties have been generally inter- 
preted in terms of a superconducting state with non-zero 
superfluid density but dominated by strong phase fluc- 
tuations of the order parameter. Such a state should 
have interesting nonlinear response characterized by per- 
sistence of superconducting nonlinearities above Tc, as 
reported before with scalar nonlinearity measurements^. 
The present study extends these results by measur- 
ing the complex harmonic voltage developed by under- 
doped cuprate superconductors when temperature is var- 
ied through the transition temperature. 

Traditionally, the microwave nonlinear response of su- 
perconductors has been investigated by using resonator 
techniques where the superconducting sample is subject 
to high microwave magnetic fields, thus making the non- 
linear effects measurable^i^. The experimental data have 
been interpreted by using various time-dependent ver- 
sions of the Ginzbug-Landau theory to estimate the or- 
der parameter relaxation time in the superconducting 
stated. After the discovery of high-Tc superconductors, 



the prospect of using these materials in microwave fil- 
ters for the wireless industry has renewed interest in the 
microwave nonlinear response. In the more recent treat- 
ments, the nonlinear effects are introduced as corrections 
to the complex conductivity and are evaluated by using a 
microscopic approach in the zero-frequency limitii^. The 
DC treatment is legitimate for the range of temperatures 
typical for the operation of high-Tc superconducting fil- 
ters (below Tc) where the superconducting order param- 
eter reacts almost instantaneously (compared to the pe- 
riod of the microwave excitation) to the applied field and 
the field screening is provided by the superfluid. The re- 
sulting field (or current density)-dependent conductivity 
is used as an input parameter for calculations of circuit 
elements in lumped-elemcnt approximations of the su- 
perconducting transmission lines and resonators^. Theo- 
retical studies addressing the operation of high-Tc super- 
conducting resonators have shown that for temperatures 
significantly below Tc, the dominant nonlinear mecha- 
nism in these devices has an inductive origin due to the 
suppression of superfluid density by the current (or ap- 
plied magnetic field)^. 

Resonator techniques have provided experimental sup- 
port for the Nonlinear Meissner Effect (NLME) at low 
temperatures in d-wave superconductors^ (enhanced by 
the presence of nodes of the order parameter on the Fermi 
surface as shown in the theoretical works of Xu, Yip and 
Sauls^ and Dahm and Scalapinoi) as well as close to Tc 
where the superfluid density is very sensitive to exter- 
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nal perturbations^!^. Despite their success, the resonator 
techniques measure parts of the sample that are often less 
than ideal, such as patterned edges, or natural edges and 
corners of single crystals. Such experiments do not pro- 
vide information about the local properties of the sam- 
ples and also usually do not provide phase information of 
the nonlinear response. This issue is relevant especially 
for high-Tc materials whose properties may vary on very 
short length scales due to their short coherence lengths. 

To overcome this limitation, a nonresonant local 
near-field microwave technique has been created to 
make spatially-resolved studies of nonlinear response of 
superconductors^iiSiii. This experimental approach is 
highly sensitive to nonlinear effects close to Tc and 
provides a high spatial resolution dictated by the geo- 
metrical dimensions of the sensing element, as demon- 
strated in harmonic measurements above an artificially- 
created grain boundar y^^i^^ . The basic idea is to ex- 
cite a highly localized current distribution at frequency 
/ on the surface of a homogeneous unpatterned super- 
conductor. Due to the nonlinear electrodynamic pro- 
cesses, harmonic (2/, 3/, etc.) signals are created 
in the material and collected by the sensing element. 
The harmonic data measured at Tc on YBa2Cu307_5 
(YBCO) thin films has been interpreted in the frame- 
work of a Ginzburg-Landau-type model where the non- 
linear source is the magnetic field(current)-dependent su- 
perfluid density n5(r, J)^. The current-suppressed su- 
perfluid density leads to an enhancement of the penetra- 
tion depth A(r, J), and consequently of the kinetic induc- 
tance, which in turn, leads to odd higher-order harmonics 
of inductive origin. The proposed model describes accu- 
rately the measured data (temperature-dependent third- 
order harmonic scalar power P3^(T)) in optimally-doped 
samples. This approach also largely avoids issues of non- 
locality that are exacerbated by current build-up at pat- 
terned edgea^iS. 

However, in underdoped samples the current- 
dependent superfluid density ns (T, J) acting alone as 
an inductive nonlinear source cannot explain the ob- 
served harmonic data. More specifically, the measured 
harmonic response P3y(T) does not turn off at tempera- 
tures above the independently determined Tc as expected 
from the model, but exhibits a tail extending signifi- 
cantly above Tc in the pseudogap regime^. It is this 
high-temperature behavior of P3/(T), including its ori- 
gin and doping dependence, that prompted the present 
study. In addition, many theoretical predictions of inter- 
esting electrodynamic properties of the pseudogap exist 
in the literatur o^'^i^^'^^i^^'^^'^^i^^ , and should be investi- 
gated. By employing a novel experimental technique, 
the harmonic response of cuprate thin films has been in- 
vestigated at temperatures close to Tc, where not only 
the magnitude of the complex harmonic voltage/power 
is measured, as in previous work, but also its harmonic 
phase. This new experimental capability also motivates a 
new field-based, finite-frequency model to explain novel 
features observed mainly in the harmonic phase exper- 



imental data. The model includes the microwave skin 
depth screening in the electrodynamics of the supercon- 
ducting state close to Tc, where the superfluid density is 
suppressed and consequently its field screening is com- 
promised. 

The paper is organized as follows: section HIl describes 
the new experimental set-up and the samples used in this 
study. Emphasis is placed on presenting in detail the 
acquisition and data processing methodology. Section 
mil presents a new theoretical model aimed at evaluat- 
ing the complex harmonic voltage measured with the ex- 
perimental set-up presented in the previous section. By 
using general electromagnetic theory and making min- 
imal assumptions about the nature of nonlinear effects, 
the model provides predictions in the temperature regime 
where ai <C cr2 and cti « (T2 ■ These limiting cases are dis- 
cussed in section HVl in conjunction with our experimen- 
tal data. The model reproduces in a semi-quantitative 
fashion some features observed in the experimental data 
as temperature is increased toward Tc. The harmonic 
data acquired on samples with various doping levels are 
discussed also in section IIVI The Discussion and Con- 
clusions sections describe the main features of the data 
and the predictions of the model. Some deficiencies of 
the model are pointed out and suggestions are made for 
future work. 



II. EXPERIMENTAL SET-UP AND SAMPLES 

The objective of the experiment (shown in Fig. [T|) is 
to locally stimulate a homogeneous superconducting thin 
film with microwave currents and measure the resulting 
nonlinear response. The microwave excitation is pro- 
vided by the internal source of a vector network ana- 
lyzer, VNA, (Agilent model number E8364B) on port 1 
at a fixed frequency / ~ 6.5 GHz in the continuous-wave 
mode, low pass filtered (to eliminate harmonics) and cou- 
pled to the sample by means of a magnetic loop probe. 
The probe is built by using commercially-available coax- 
ial cable (UT034) where the inner conductor has been 
soldered to the outer one^iiiSiii. This results in a semi- 
circular loop with inner radius of roughly I65/im and an 
outer one of 365/xm (see inset of Fig. [T]). The loop has 
been further mechanically polished at the outer radius in 
order to bring the microwave current, flowing in a thick- 
ness dictated by the skin depth at the inner radius of the 
loop, closer to the sample (~ 100/xm), thus improving the 
loop-to-sample coupling and allowing the operation of the 
apparatus at lower input power levels. The signal origi- 
nating from the sample comes back through the probe, is 
high-pass filtered to suppress the microwave power at the 
fundamental frequency / and examined with the VNA 
in the frequency-offset mode (VNA-FOM) by tuning the 
receiver on port 2 in a narrow frequency range (I Hz) 
centered on the harmonic of interest (3/ in this case). 

The sample is placed in a cryogenic environment whose 
temperature is controlled between 78 and 100 K with an 
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accuracy of 0.1 K and for each temperature a trace is ac- 
quired from the VNA-FOM and stored on a computer for 
further analysis. The samples are unpatterned homoge- 
neous c-axis oriented YBa2Cu307_5 (YBCO) thin films 
deposited on SrTiOs or NdGaOs substrates by pulsed 
laser deposition (PLD). The oxygen content of the sam- 
ples has been adjusted by annealing in various oxygen 
pressures and at different temperatures in the PLD cham- 
ber resulting in critical temperatures T^'^ in the range 
of 52 to 90 K, as evaluated from AC susceptibility mea- 
surements (performed at 120 kHz). Despite their small 
thickness (~ 50 nm) the sample superconducting qual- 
ity is very good as revealed by the narrow peaks of the 
temperature-dependent imaginary part of the magnetic 
susceptibility whose full width at half maximum ST^'~^ 
are given in Table HI 



Sample 


T^^ [K] 


5r,^^ [K] 


7~S 


ATM.m [K] 


Substrate 


SI 


88.9 


0.3 


6.84 


0.2 


NdGaOa 


S2 


86.6 


1.0 


6.82 


0.5 


NdGaOa 


S3 


74.0 


0.9 


6.76 


0.4 


NdGaOa 


S4 


62.0 


0.55 


6.69 


0.6 


SrTiOs 


S5 


52.0 


1.1 


6.53 


1.0 


SrTiOa 



TABLE I: Sample properties: critical temperature and 
transition width ST^'~' determined from AC susceptibility 
measurements, the doping level 7 — 6 estimated from T^'~^ , 
the difference between the temperatures where the extreme 
values of the harmonic phase and magnitude occur, ATM,m, 
and the sample substrate. 

After a power calibration, performed according to the 
manufacturer's instructions^^, the VNA-FOM measures 
the absolute power level at frequency 3/, -P3/, incident 
on port 2 similar to a spectrum analyzer (" spectrum an- 
alyzer mode" ) . 

To perform phase-sensitive detection of the 3/ har- 
monic voltage incident on port 2, the VNA-FOM requires 
a reference signal at the same frequency as the signal to 
be analyzed. For the measurements reported here, the 
reference signal is provided by an additional microwave 
circuit, called the reference path (see Fig.[T|) that converts 
some microwave power at frequency / generated by the 
internal VNA source into microwave power at frequency 
3/. The fundamental / is fed into a comb (harmonic) 
generator (Herotek, model number GCA 2026A-12) fol- 
lowed by two band-pass filters designed to suppress (at- 
tenuation of 80 dB) the fundamental and all harmonics 
except for 3/, thus resulting in a clean 3/ signal that 
serves as reference, U^'t'^ , for the phase-sensitive har- 
monic detection. The VNA-FOM traces represent the 
complex ratio of the voltage from the sample (whose tem- 
perature T is varied inside the cryostat) , to that from the 
reference path, U3/(T) = U^'^"'^^'^ (T) /Ul^y at the plane 
of the VNA's port 2, evaluated at the frequency points 
scanned by the VNA-FOM receiver within the 1 Hz span 
window centered on 3/ ("vector signal analyzer mode"). 
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FIG. 1: (Color online) The experimental set-up for the phase- 
sensitive harmonic measurements. A signal at frequency / 
is generated by the internal VNA source, low-pass filtered 
and delivered to the sample inside a cryogenic environment. 
The microwave signal reflected by the sample and containing 
harmonics of the incident power is high-pass filtered to re- 
move the fundamental / and measured on VNA port 2. For 
phase-sensitive detection a reference circuit, converts some 
microwave power at frequency / into power at 3/ by means 
of a harmonic generator and two band-pass filters, and feeds 
the resulting 3/ signal back into the VNA as reference signal. 
Inset: drawing of the magnetic loop probe above the sample. 



Since during an experiment the reference path and most 
of the microwave circuit are at room temperature (only 
about 10 cm of coaxial cable is inside the cryostat, how- 
ever not in physical contact with the cold plate), it is le- 
gitimate to attribute the temperature dependence of the 
measured relative harmonic voltage u^f{T) entirely to 
the temperature-dependent nonlinear effects in the sam- 
ple. Since the microwave circuit is operated at low mi- 
crowave power (~ to -1-9 dBm at VNA port 1), the back- 
ground noise at 3/ is dominated by the intrinsic noise of 
the VNA (~ —140 dBm when P^f is measured without 
an input on VNA port 2). This is shown in Fig. ^ bot- 
tom plot, where the nonlinear signal from the sample are 
at the noisefloor (T=83.4 and 91.5 K). 

The VNA-FOM is not designed for absolute phase 
harmonic measurements, but for relative ones since 
the VNA-FOM only indicates how the harmonic phase 
changes as the sample properties change from one tem- 
perature to another, provided that the reference sig- 
nal U^f'^ is stable during the measurements. Conse- 
quently, the temperature-dependent VNA-FOM phase 
data 03/ (T) = ci>-™f'-(r) - - ^offset, are offset 
by a temperature-independent unknown amount {^^Y + 
^offset) originating from the phase winding in the com- 
ponents of the microwave circuit (coaxial cables and fil- 
ters) and the phase relationship between the fundamen- 
tal and the harmonics generated by the comb generator. 
The phase shift + ^of fset is evaluated by using the 
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FIG. 2: (Color online) Examples of VNA-FOM traces ac- 
quired on a YBCO thin film (Sf) in a frequency range cen- 
tered on 3/ — f 9.47 GHz. Top and middle plot: the phase and 
magnitude of U3f{T) acquired in a phase-sensitive measure- 
ment; bottom plot: absolute harmonic power data P3/(r), 
acquired in "spectrum analyzer mode". 

predictions of the theoretical model presented in Section 
mil for the limiting case of low temperatures (T ^ Tc) 
where the harmonic phase is 7r/2, with a minimum of 
assumptions and in agreement with other experimental 
observations^. 

To increase the signal-to-noise ratio of the VNA-FOM 
traces, 8-10 averages were performed on the VNA be- 
fore transferring the trace to the computer. Such av- 
eraged traces acquired at three representative temper- 
atures (below, above and around T^'^) are shown in 
Fig. [2| the complex phase and magnitude of the har- 
monic voltage (upper and middle plots) and the abso- 
lute harmonic power levels (bottom plot). The 1 Hz 
span P^f VNA-FOM traces look similar to small-span 
traces acquired with a spectrum analyzer: at T=89.6 
K (close to T^'-^) the temperature-dependent harmonic 
power P^f{T) reaches its maximum in agreement with 
other authorsMiSiio, while at T=83.4 K and T=91.5 K 
(below and above T^^ respectively) the traces are flat 
and at the noisefloor (—135 to —140 dBm). Complex 
M3/(r) VNA-FOM averaged traces acquired at T=83.4 
K and T=91.5 K, exhibit a large scatter both in phase 
and magnitude, while |M3j^(T)| reaches its maximum at 
T=89.6 K, in agreement with the P3/ traces. 

To extract the temperature dependence of the relative 
harmonic voltage u^f{T) from the VNA-FOM traces, 1 
to 3 central points from each trace (see Fig. [2]) are aver- 
aged, and the resulting complex magnitude and phase are 
plotted vs. temperature in Fig. [3l To quantify the data 
spread in a U3f{T) VNA-FOM trace and the reliabihty 
of averaging the complex traces, the standard deviation 
of the phase data STD^^^ is evaluated from each of the 
averaged 1 Hz span traces and represented together with 
the phase data in Fig. [31 The temperature-dependent 
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FIG. 3: (Color online) Phase-sensitive harmonic data ac- 
quired on a YBCO thin film (SI), a): temperature-dependent 
magnitude, |M3/(r)|. b): temperature-dependent phase 
(03/ (T) green solid line) and the standard deviation of the 
11-point traces acquired at each temperature (blue dashed 
line). Magenta arrows show: T=83.4, 89.6 and 91.5 K. 



STD^^j can be used to select a temperature range where 
03 / (T) data can be considered reliable by imposing that 
STD^^^ does not exceed a certain threshold. The tem- 
peratures associated with the traces from Fig. [5] are indi- 
cated in Fig. [3] with magenta arrows together with T^'~^ 
and ST^'-^ from AC susceptibility measurements. 

The samples have been measured by using various in- 
put frequencies (6.45 to 6.55 GHz) and power levels (0 to 
+9 dBm) and with the microwave probe placed at several 
locations above the samples, all with consistent results. 
For the range of microwave input power levels employed 
in this work, the harmonic data suggest that the mi- 
crowave probe does not induce a significant amount of 
heating in the sample surface, which would be indicated 
by a shift of the maximum of |{(3j^(T)| and of the mini- 
mum of (/)3 / (T) respectively to lower temperatures when 
the microwave power is increased. 

Qualitatively, the magnitude of the harmonic voltage, 
|u3/(T)|, reaches a maximum at a temperature Tm close 
to T^'~^ in agreement with results from the literature 
where P3f{T) is reported to reach a maximum^i^'^'i^. 
One novelty of the results presented here is that the 
complex phase decreases smoothly and reaches a min- 
imum at Tm as the temperature is increased. For the 
samples investigated here a consistent trend has been 
observed: for near-optimally doped samples the temper- 
atures associated with the two extrema of magnitude 
and phase almost coincide, while for underdoped sam- 
ples the harmonic phase tends to reach its minimum at 
higher temperatures (Tm > Tm)- To quantify this trend, 
ATm.™ = Tm — Tm has been evaluated for all samples 
from measurements at various input frequencies and mi- 
crowave power levels and is given in Table HI There is a 
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general trend of increasing ATM,m with increased under- 
doping. 



III. MODEL 

The recent theoretical treatments of microwave non- 
linear effects in superconductors are mostly restricted 
to resonant configurations and their equivalent lumped- 
element circuit approximations^. Due to the non- 
resonant nature of the near-field microwave microscope 
employed in this work and since the sample is not part 
of a transmission line structure, a field-based analyti- 
cal model is more appropriate to capture the essential 
physics. Such an approach has been proposed for tem- 
peratures below Tc, where the authors considered only 
the nonlinear effects caused by the field dependence of 
the imaginary part of the complex conductivityi^. 

The analytical model presented in this section is an ex- 
tension of that proposed by Pestov and co-workers^*' at 
finite frequencies in an attempt to explain the temper- 
ature dependence of the harmonic phase data acquired 
with the near-field microwave microscope described in 
the previous section. The goal of the model is to pro- 
vide an expression for the complex-valued temperature- 
dependent harmonic voltage U^Y^^^'^ {T) induced in the 
near-field antenna by the microwave screening current 
distribution from the nonlinear superconducting sample. 
This is achieved in three steps: first, the magnetic vector 
potential created by the excitation current in the mi- 
crowave probe, is calculated at the sample surface, 
then by using a nonlinear generalization of the constitu- 
tive London equation for superconductors, the harmonic 
content of the screening current induced in the sample 
is evaluated, and in the third step the harmonic voltage 
induced in the near-field probe is found. 

In order to preserve a higher level of generality, the 
nonlinear effects are introduced here as phenomenological 
corrections to both the real and imaginary parts of the 
low-power, linear-response complex conductivity of the 
sample a = ui — ia2, where ai^2 are positive definite 
(cti,2 > 0): 

ai,2(r, A/) = (71,2 ± ^ + • • j , A/ « Aia (1) 

Af is the vector potential associated with the microwave 
excitation at the fundamental frequency / and the 
temperature-dependent nonlinear vector potential scales 
Ai^2 quantify the nonlinear effects in the real and imag- 
inary components of the complex conductivity, respec- 
tively. We assume local electrodynamics in this model. 
Note that Ai^2 can model a wide variety of nonlinear 
sources, including the nonlinear Meissner effect, vortex 
motion, weak links, etc. These corrections are valid when 
Af <^ ^1,2 J similar to other phenomenological descrip- 
tions of nonlinear effects in the literature^-ii^iSi . This 
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FIG. 4: (Color online) The geometry of the analytical prob- 
lem. The near-field microwave antenna is approximated by a 
filamentary current wire extending infinitely in the z direction 
aX X = a, y = Q. The sample with thickness do extends in- 
finitely in the zOy plane and has its bottom surface at a; = 0. 



treatment also implicitly assumes that the complex con- 
ductivity of the superconductor reacts instantaneously to 
changes of the probing vector potential Af. At temper- 
atures very close to Tc, where the condition ^4/ <C Ai,2 
might be violated and higher order terms should be in- 
cluded in the expansion ([1]), or when order parameter 
relaxation times become comparable to the microwave 
period, the present formalism might not be applicable. 
In a qualitative picture the superfluid density ns is sup- 
pressed by the microwave excitation ("-" sign in Eq.([T|) 
for (72 ) and "converted" into normal fluid ("-I-" sign in 
Eq.(P) for (Ti). 

For the analytical treatment of the problem, the near- 
field probe is approximated by a filamentary current wire 
parallel to the z axis located at y = and x = a whose 
current density is modeled by the 5 Dirac distribution: 
jY^i^jV) = If^{^ ~ o,)^{y)i while the sample extends 
infinitely in the zOy plane of the Cartesian frame with 
the bottom surface at x = (see Fig. |4]). Since the 
wire-to-sample separation, a, is much smaller than the 
wavelength of the microwave excitation, the magnetic 
vector potential Af outside the superconducting sam- 
ple satisfies the equations of magnetostatics: above the 
sample V^A/ — — /ioj/^*(2:, 2/) while below the sample 
V^A/ =0. Inside the superconducting sample the spatial 
variation of the magnetic vector potential A f and current 
density in the film jf um in the x direction is governed by 
the length scale associated with the inverse of the wave 
vector 17""'^ I (7^ = A^'^ -I- 2i5^^{f) which at low temper- 
ature T ^ Tc is governed by the penetration depth A 
while in the normal state by the skin depth evaluated at 
the probing frequency /, 5sk[f))- For details on the su- 
perconducting screening at finite frequencies, see for ex- 
ample the work of Clem and Coffey22i. Since the sample 
thickness is smaller than both the zero-temperature pene- 
tration depth as well as the microwave skin depth within 
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the investigated temperature range {do <C Aq, 
Af and jfum are assumed uniform within the thickness 
do. At microwave frequencies the displacement current 
is neghgible with respect to the conduction one, there- 
fore Ampere's law in integral form for the closed loop L 
(shown in Fig. |4|) reads: 

Bf{x = do + 0) - Bf{x - 0) = fiojfamdo (2) 

where jfumdo is the sheet current. In the limit do — > of 
sample thickness much smaller than all length scales in- 
volved in the problem, and taking into account that the 
current density inside the film obeys London's law for 
finite-frequency: jfUm = — 7^^///^o (here we assume lo- 
cal electrodynamics), the equation for the magnetic vec- 
tor potential, can be written in a closed-form for the en- 
tire spacei^: 

-W^Af{x,y) + X;j\Af{x,y)S{x) = fioIf6{x~a)6{y) (3) 

where Ae// — l/((io7^) represents a generalized finite 
frequency effective penetration depth. 

To integrate Eq.®, the nonlinear effects in cri^2 are 
neglected in this step {Ai^2 ~* oo) and the equation is 
Fourier-transformed. After solving for Af{kx, ky) and in- 
tegrating with respect to kx, Af{ky) readsi^: 



Xeff exp(-|fc^|( 

l + 2Xeff\ky\ 



(4) 



For the experimental setup described here, the sample- 
to- wire separation, a is determined by the diameter of the 
inner conductor of the coaxial cable (a 100/im). In the 
long-wavelength approximation, a exceeds both length 
scales contained in Ag// (A and Sgkif)) and consequently 
the term 2Xeff\ky \ can be neglected in the denominator 
of Eq. ^ , allowing a closed- form expression for the vector 
potential generated by the current wire: 
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The finite-frequency nonlinear generalization of the lo- 
cal London constitutive relationship jfUm — js + jn — 
((Ti — ia2)E — —Lo{i<Ti + '^2)Af (where the electric field 
E = —dAf/dt and Aj ^ exp{+iujt) with cu = 27r/) is 
obtained by replacing the linear-response complex con- 
ductivity (Ti^2 with its phcnomenological nonlinear ex- 
pressions from Eq.(IT]): 



jfUrn ^ -UJ<J2\l - -MAf 



iw^ifl + ^W (6) 



This approximation is valid under the limited condi- 
tion Af <^ Ai 2 and shows that the current density con- 
tains a component at frequency / and another component 
at frequency 3/ which represents the source of the mea- 
sured harmonic voltage at frequency 3/. The nonlinear 
component at frequency 3/ in the total current density 
jfiim, jaf, is separated from the A^j terms by consider- 
ing the time dependence Af cos{ujt) and using the 
trigonometric relation cos^ ujt — {cos 3u!t + 3 cos ujt) / A: 
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The current distribution flowing in the sample and 
having a harmonic 3/ time variation generates a vec- 
tor potential A^ f in the entire space and induces a volt- 
age in the near-field probe. In order to evaluate the 
induced voltage at 3/, U^J^^^'' , one has to calculate 
the vector potential at the location of the wire. This 
is accomplished by using the equivalence principle from 
electromagnetism^ where a current with frequency 3/ 
flowing through the wire jg^' = l3fS{x — a)S{y) generates 
the magnetic vector potential on the sample surface given 
by Eq.® with the appropriate substitution / 3/. 
Equivalently, a current distribution j^f in the sample 
given by Eq.® generates a vector potential A3^(x,y, z) 
in the entire space. The equivalence principleSS, reads: 



dVjifix, y, z)A^f(x, y, z)^ dVj3f(x, y, z)A3f{x, y, z) (9) 



with the integrals evaluated over the entire space. By 
using the filtering properties of the Dirac delta function 
and since all the z = constant planes contain the same 
field and current configuration due to the symmetry of 
the problem, the vector potential at the location of the 
wire reads: 
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X-^ + 2i5-^'{f)J X-^ + 2z5;,'{3f) 



(10) 



where 6sk{f) and (5sfc(3/) represent the skin depth evalu- 
ated at frequency / and 3/ respectively. These two quan- 
tities differ by a factor of V3 in the ordinary skin effect 
regime, and in order to simplify the calculations, the fol- 
lowing approximation will be used: (5sfc(3/) ~ <^sfe(/) = 

Ssk- 

The electric field induced in the wire at frequency 3/, 
E^f = —dA^^{a,0)/dt, is used to evaluate the voltage 
induced in a probe of length Iq: 



or can be expressed in terms of the linear-response length 
scales A and Sgk- 
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In terms of conductivities, the induced vohage reads: 
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(12) 



Equations pT|) and (fT2|) have been deduced in an an- 
alytical, field-based approach, as opposed to most recent 
models in the literature which use lumped-element de- 
scriptions for the superconducting devices operating at 
microwave frequencies. Several features can be noted: 
the harmonic voltage magnitude scales with the exci- 
tation current as \U^'^™'^'''^\ ~ /j, the nonlinear effects 

are easier to measure in thin films {\U^j™'^'''^\ ~ ^cT^) 
and for small antenna-to-sample geometric separation 
(jC/g^™''''^! ~ a~^), all in agreement with experimental 
data in the literature, as well as models. 

The model provides an estimate for the complex- valued 
harmonic voltage induced in the near-field probe by the 
screening current flowing on the sample surface. Since 
the VNA-FOM measures the harmonic voltage from the 
sample with respect to that from the reference path at 
the plane of VNA's port 2, the measured data must be 
phase-shifted by an amount + $0/ /set (see Section 
HI)) , which is equivalent to moving the measurement plane 
from VNA's port 2 to the near-field antenna. Such a 
translation is effectively accomplished by examining the 
limiting case T ^ Tc of Eq. to evaluate the required 
amount of phase shift. 



IV. 



ANALYSIS 



A. Harmonic Phase vs. Temperature 

Before comparing the results of the mathematical 
model with the experimental data, it is useful to exam- 
ine the case of nonlinearity at low temperatures, which 
has also been investigated by other authors in a reso- 
nant configuration^!. For T ^ Tc, when the contribu- 
tion of the normal fluid to the electrodynamics of the su- 
perconducting state is small, the induced voltage given 
by Eq. (jl2p can be expanded in a power series around 

CTl/tT2 = 0: 
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Two possible scenarios emerge from this picture de- 
pending on the ratios of conductivities and that of vec- 

tor potential scales. For the case 1 3> ^ ^ the 

harmonic voltage ?7|?™^''^(a, 0) has a zero phase and de- 
pends on the a\ nonlinearity characterized by the vector 



potential scale A\: 
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On the other hand, for the case — <C ^4- the first term 



in the power expansion Eq. (|13p is dominant and repre- 
sents a pure inductive- like nonlinear response. The har- 
monic voltage {7|^™P''^ q) depends only on the (T2 non- 
linearity, characterized by the nonlinear vector potential 
scale A2 , as in most of the treatments of superconductor 
nonlinear response^'^ i^^'^"'^ : 
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This result does not rely on any assumptions about the 
temperature dependence of conductivity correction fac- 
tors Ai,2 and is in good agreement with the experimental 
data of Booth et al. who reported a purely inductive har- 
monic response in YBCO thin films at 78 K^^, well be- 
low Tc « 90 K. Additionally, theoretical models describ- 
ing the operation of high-Tc microwave filters at tem- 
peratures below Tc, predict a purely inductive nonlinear 



response (consistent with our case fj' ^ 



» 1) 



however do not exclude the possibility of resistive nonlin- 
ear effects close to T(M. Equation ([T5]) shows that in the 
superconducting state where <C 1, the experimental 
setup is not sensitive to a possible ai nonlinearity; more- 
over, even if the cti nonlinearity would dominate the CT2 
one (— <C ^ <C 1), this effect would not be detectable 

^ (72 A^ ' ' 

due to the modest contribution of g\ in the field screening 
process. 

The harmonic phase data reported here exhibit an al- 
most fiat plateau at low temperatures where the magni- 
tude starts to go above the noisefloor (see, for example, 
FiglS]). This observation, together with the prediction of 
Eq. suggests that the harmonic response in that tem- 
perature range is characterized by a 7r/2 phase. By using 
this result, the relative phase data acquired by the VNA- 
FOM on all samples have been corrected over the entire 
temperature range by adding a temperature-independent 
phase offset (labeled ^^-^Y + $0/ /set in section |Tl| to en- 
force the condition (f>3f{T<^ Tc) ~ tt/2. 

As temperature is increased toward Tc, the ratio g\I 02 
increases and the in-phase component starts to become 
significant (the second term in square brackets in Eq. (jl3p 
which is real and positive) while the out-of-phase com- 
ponent (the difference of the first and third terms in 
Ea. (fT3l) ) is gradually reduced. This behavior is consis- 
tent with the data from all samples, (see, for example, 
Fig|3]) showing that the third harmonic phase angle ro- 
tates clockwise from 7r/2 as Tc is approached from below. 
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FIG. 5: (Color online) Temperature-dependent phase- 
sensitive third-order harmonic voltage data U3 / acquired on a 
YBCO (SI) thin film represented in the complex plane. The 
arrows indicate the evolution of the complex data as tempera- 
ture increases from Start(T=86.4 K) to End(T=90.4 K). Only 
the low-noise data are presented, here in arbitrary units. The 
dots represent the predicted phase at a temperature for which 
CTi = (J2 and either Aj < A2 or A^. 



Equation for tlie liarmonic voltage includes tlie ra- 
tios cri/a'2, ^1/^2 J ^ii'i ^2 whose temperature depen- 
dence must be known in order to model the experimental 
values of the harmonic voltage magnitude and phase. For 
a semi-quantitative discussion it is important to examine 
the temperature-dependent phase of the third harmonic 
voltage by considering only the functional dependence of 
the temperature-dependent terms in Ea. (|12p : 
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The behavior of the harmonic voltage in the complex 
plane is dominated by the last term in Eq. (|16p due to its 
4*'' power. For superconductors, in a mean-field approx- 
imation, the ratio of conductivities (Jijo^ is essentially 
zero at T ^ Tc and generally increases towards infinity 
as the temperature approaches Tc; therefore the temper- 
ature dependence of the complex phase associated with 
the last term in Eq. ()16p is at low temperature and exe- 
cutes a full 360 degrees clockwise rotation in the complex 
plane as temperature is increased toward T^. This pre- 
diction for the sense of rotation is in agreement with the 
experimental data up to the temperature T^ where 03/ 
goes through its minimum. 

The harmonic voltage U'^^^^^ia,^) given by Eq. fT^ 
can be expanded in a power series around a\ j 02 = 1 and 
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In the limiting case a\ — the complex harmonic volt- 
age \J^^^^^(a, 0) Hes in the 3'''' quadrant of the complex 
plane with negative real and imaginary parts. To check 
this theoretical prediction, the experimental data shown 
in FigO have been represented in the complex plane as 
Re(M3/) vs. Im(u3/) in FiglSl after offsetting the phase 
data to enforce the condition (^3/ « 7r/2 at the lowest 
temperature where the signal-to-noise ratio is acceptable. 
Also in FigOthe two extreme cases of A\ ^^2,0-1 = a^, 
and A\ :s> Al,(Ti — (J2 have been represented as two dots. 
Common to all samples from Table HI the harmonic phase 
data exhibits a non-monotonic behavior, decreasing from 
7r/2 at low temperatures, reaching a minimum inside the 
4*'' quadrant at a temperature T^, and increasing back 
toward 7r/2. The experimental data acquired with all the 
samples from TableUdo not reach the 3'"'' quadrant of the 
complex plane, as Eq. pT)) predicts. Note that the model 
does not take into account the finite order parameter 
relaxation time which, close to the critical temperature 
(i.e. in the regime ui « (72), becomes comparable to the 
microwave probing period. 

Experimental data acquired at low frequency (1 kHz) 
by Mawatari and co-workers^^, although in a different ex- 
perimental configuration that enhances the electromag- 
netic response associated with vortex motion, exhibits 
a similar trend. The 3'''' order harmonic phase goes 
through a minimum of roughly — 7r/2 when superconduc- 
tivity is gradually weakened by the application of an ex- 
ternal magnetic field, whereas for the data presented here 
the superconductivity is suppressed by increasing tem- 
perature. The microwave harmonic measurements pre- 
sented here have been performed on a time scale roughly 
7 orders of magnitude shorter than that employed in the 
work of Mawatari et ali^, thus they are more prone to 
explore the regime where the superconducting order pa- 
rameter cannot oscillate in phase with the probing field. 
This could prevent the harmonic phase data from reach- 
ing the 3'''^ quadrant of the complex plane. The mini- 
mum of the harmonic phase observed in the microwave 
harmonic data could indicate the onset of the regime 
where the dynamics of the superconducting order pa- 
rameter becomes slower than the period of the probing 
electromagnetic field. To account for this regime and 
to explain third-order harmonic data acquired in gapless 
superconductors, Amato and McLean^ proposed a cor- 
rection term (1 — 2iu!Tji)~^ to the microwave harmonic 
field Bjf, where Tf({T/Tc) is the order parameter relax- 
ation time that diverges at Tc. Such a correction term 
gives rise to a counter-clockwise rotation of the harmonic 
voltage in the complex plane, in agreement with our ex- 
perimental data at temperatures above T^. 
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FIG. 6: (Color online) Comparison of the third-order har- 
monic voltage magnitude data versus temperature acquired 
on the samples listed in Table |I] 



B. Harmonic Signal vs. Doping 

A comparison of the magnitude and phase data ac- 
quired on all samples from Table H] is shown in Figl6]and 
[21 Only the relatively low-noise data have been presented 
in these figures: for magnitude data, the noisefloor is 
given by the sensitivity of the VNA-FOM, while for phase 
data only the temperature ranges where the phase stan- 
dard deviation STDg,.^^ is less than 0.2 radians have been 
selected. 

The examination of both Fig. [6] and Fig. [7] reveals that 
the maximum of the harmonic voltage magnitude occurs 
at a temperature Tm, lower than that associated with 
the minimum of the complex phase, T^- The differ- 
ence between these two temperatures, denoted ATjv/.m — 
Tm — Tm, follows a consistent trend as indicated in Ta- 
ble H] and shown in the inset of Fig. [7] (i.e. increasing 
in the more underdoped samples), and is not correlated 
with the broadening of the superconducting transition in 
underdoped samples due to the annealing process (see 
Table [T] for ST^'~'). This feature has been observed in all 
measurements despite the different input power levels, 
probing frequency or microwave antenna location above 
the sample. Unfortunately, due to the unavailability of 
ai^2 and A1.2 theoretical temperature- and doping depen- 
dences, the observed trend of ATM,m with doping could 
not be accounted for by the model. However, future theo- 
retical models could be checked against this experimental 
observation. 

The harmonic voltage magnitude data from Figl6]show 
that in near optimally-doped samples, the harmonic volt- 
age drops sharply into the noisefloor at a temperature 
slightly above Tm, while in underdoped samples, the |m3/| 
peak occurs at a temperature below Tm and the \u3f\ 
temperature dependence extends more above Tm- This 
general trend is in agreement with the observation of Lee 
et at, who examined the microwave harmonic power P^f 
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FIG. 7: (Color online) Comparison of the phase of the third- 
order harmonic voltage data versus temperature acquired on 
the samples listed in Table HI Inset: doping dependence of 

ATM.m. 

reflected by YBCO thin films with various doping levels^. 
An interesting feature can be noticed in the third-order 
harmonic voltage phase data (Fig.[7|): in optimally-doped 
samples when the magnitude |u3/(T)| reaches the noise- 
floor the harmonic phase recovers its value from the su- 
perconducting state {+Tr/2) after exhibiting a minimum 
at Tm, whereas in the most underdoped samples, the re- 
covery is not complete. 

Some features in FiglHland Fig[7]can be explained qual- 
itatively by using the argument of Amato and McLean^ 
regarding the regime where the superconducting order 
parameter cannot react instantaneously to the microwave 
excitation. The weakening of the harmonic phase reversal 
and the broader extent of the harmonic magnitude ob- 
served in our harmonic data measured in the underdoped 
samples points to a slower divergence of the order param- 
eter relaxation time in underdoped samples compared to 
their optimally-doped counterparts. Such a scenario is 
compatible with stronger and longer-lived superconduct- 
ing fluctuations in underdoped cuprates. The extension 
of the nonlinear response above Tc is consistent with An- 
derson's picture in which the pseudogap phase has the 
electrodynamic properties of a superconductor, but with 
a current-current correlation function that decays with 
a finite time r, and a diamagnetic susceptibility that is 
nonlinear—. 



V. DISCUSSION 

Overall, the model offers a unified picture of microwave 
nonlinear effects originating from both the real and the 
imaginary parts of the conductivity. The model repro- 
duces semi-quantitatively the trends observed in both the 
magnitude and phase of the harmonic voltage acquired 
with the phase-sensitive nonlinear near-field setup. The 
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lack of theoretical models for the temperature- and dop- 
ing dependence of ai 2 and the divergence of the order 
parameter relaxation time close to the critical tempera- 
ture hinders a more detailed comparison of data with the 
predictions of the model and the extraction of the Ai^2 
temperature- and doping dependence. 

The analytical treatment points out the significance of 
the ai nonlinearity in addition to that associated with 
(72. Most of the recent treatments of the nonlinear ef- 
fects in cuprates close to the critical temperature Tc ig- 
nore the influence of the normal fluid and consequently 
assign the observed nonlinear effects entirely to the (T2 
nonlinearity^iiS. Well below T^, the observed nonlinear 
behavior is due to the superfluid, as shown in the theoret- 
ical works of Dahm and Scalapinc^, confirmed by the ex- 
perimental work of Bootb^ and also shown in our model 
for the case of low temperatures (cri ^ (12). However, 
the microscopic model of Dahm and Scalapino does not 
exclude the possibility that in close proximity to Tc non- 
linear mechanisms due to the normal fluid might become 
important^*. 

In the model, the relaxation time of the order param- 
eter was assumed much shorter than the microwave pe- 
riod, i.e. the order parameter reacts instantaneously with 
changes in the external probing field. This assumption is 
valid only up to temperatures very close to T^^. The new 
data presented here correlated with those of Mawatari et 
ai2£ acquired at much lower excitation frequency suggest 
that neglecting the dynamics of the superconducting or- 
der parameter restricts our analysis to temperatures be- 
low Tm where the phase of the harmonic voltage reaches 
its minimum in the 4*^* quadrant of the complex plane. 

The major benefit of the experimental technique comes 
from the localized nature of the microwave excitation 
and the ability to measure the temperature-dependent 
complex harmonic voltage. Therefore, the nonlinear mi- 
crowave response can be investigated in as-grown super- 
conducting samples, free from potential defects caused by 
patterning. Additionally, the harmonic response can be 
measured at various locations in a homogeneous sample, 
thus ensuring that the response does not depend on loca- 
tion or some peculiar feature of the sample, such as edges, 
corners, grain boundaries or defects of fabrication^^. 

The ability to measure both the magnitude and the 
phase of the harmonic voltage allows a more complete de- 
scription of nonlinear effects as a function of doping level. 
The experimental data show that the maximum of har- 
monic voltage magnitude and the minimum of phase oc- 
cur at slightly different temperatures, Tm and Tm respec- 
tively. The doping dependence of the nonlinear response 
has been quantified by defining ATpi.-m — Tm — Tm 
and monitoring its variation with doping level 1 — 5. 
In almost optimally-doped samples the harmonic voltage 
magnitude exhibits a sharp maximum very close to the 
temperature where the harmonic phase reaches its mini- 
mum, Tm, then drops abruptly to the noisefloor. In the 
more underdoped samples, the maximum of the harmonic 
voltage magnitude occurs at a temperature lower than 



Tm, thus in the superconducting state, but the harmonic 
response extends above Tm into the pseudogap phase. 
This doping-dependent nonlinear response could be due 
to enhanced Cooper pair lifetimei^ii^, to superconduct- 
ing fluctuations in the pseudogap phase of underdoped 
cupratesi^ii^iiiii^ or perhaps to a transition from a pure 
d-wave order parameter to a, d + s order parameter—. 
Another source of nonlinear response above Tc could be 
vortex-like excitations in the pseudogap phase^i. 

By a proper choice of the complex conductivity ui 2 
and its nonlinear corrections ^1,2, one could incorporate 
microscopic details (symmetry of the order parameter, 
shape of the Fermi surface, effects due to the quasiparti- 
cles at the nodes of the gap, etc.), anisotropy effects of the 
in-plane conductivity and various possible sources of non- 
linear behavior, such as vortex motion, weak links and 
defects due to sample fabrication and annealing. Further 
investigations of nonlocal electrodynamics^ should also 
be explored along with more microscopic models. 



VI. CONCLUSIONS 

A new experimental technique is presented where a 
vector network analyzer in frequency offset mode is used 
to acquire the harmonic nonlinear response of homoge- 
neous thin superconducting films to microwave current 
excitation. The phase-sensitive harmonic detection tech- 
nique provides an additional piece of information, com- 
pared to previous investigations of nonlinear effects in 
the superconducting state: the phase of the harmonic 
voltage at temperatures close to Tc. The third-order 
harmonic phase gradually decreases from 7r/2 in the su- 
perconducting state, reaches a minimum close to Tc and 
recovers back to 7r/2, at least for samples near optimal 
doping. In the underdoped samples the phase does not 
recover completely to 7r/2 in the normal state, but to 
roughly 0.5 radians. In all samples used in this study 
the harmonic magnitude exhibits a maximum, as ob- 
served by other investigators^iiS and the magnitude max- 
imum occurs at a temperature Tm below that associ- 
ated with the minimum of phase, Tm- A consistent trend 
with oxygen doping has been found, where the difference 
ATm, ~Tm increases in the more oxygen-deficient 

samples. 

An analytical finite-frequency field-based model of the 
nonlinear microwave response of superconducting thin 
films in a near-field microwave experimental configura- 
tion is presented. The interplay of inductive and resistive 
nonlinear effects arises naturally in the model, being an 
improvement with respect to previous models from the 
literature, which treat the two types of nonlinear behav- 
ior separately. The description is field-based as opposed 
to lumped-element-based and it introduces the nonlin- 
ear effects in a phenomenological fashion as deviations 
of conductivity from its linear-response value. The new 
model is in agreement with experimental data and other 
models from the literature in the limiting case of low 
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temperature where the field screening is due to the su- 
perfluid, and reproduces some key features observed in 
the data acquired with the new apparatus. More specif- 
ically, at low temperatures it predicts a harmonic phase 
of 7r/2 with a nonlinear behavior originating mainly from 
the (72 nonlinearity. The model also shows, in agree- 
ment with our experimental data, that as temperature 
increases toward Tc an in-phase, resistive-like, compo- 
nent becomes significant, thus "rotating" the third-order 
harmonic voltage clockwise in the complex plane. In the 
regime where the superconducting order parameter re- 
laxation time diverges, the counter-clockwise rotation of 



the harmonic voltage in the complex plane can be quali- 
tatively reproduced by a phenomenological modification 
of the model. 
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